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Abstract 

Current data implies three simple forms of the neutrino mass matrix, each corre- 
sponding to the conservation of a non-standard lepton charge. While models based 
on Lg and Lg — — Lr are well-known, little attention has been paid to — L,-. A 
low energy mass matrix conserving — Lj- implies quasi-degenerate light neutrinos. 
Moreover, it is /x-r symmetric and therefore (in contrast to Lg and Le — — Lr) 
automatically predicts maximal atmospheric neutrino mixing and zero C/gS- A see- 
saw model based on — is investigated and testable predictions for the neutrino 
mixing observables are given. Renormalization group running below and in bet- 
ween the see-saw scales is taken into account in our analysis, both numerically and 
analytically. 
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1 Introduction 



Impressive steps towards identifying the structure of the neutrino mass matrix rrii, have 
been taken in recent years p. In the charged lepton basis rrii, is diagonalized by the PMNS 
matrix j2] U defined via 

f/'^ mj, [/ = diag(mi, m2, ma) . (1) 

It can be parametrized as 

/ C12C13 S12C13 sise"*'' \ 

U = I — C23 S12 — S23 Si3 C12 e*^ C23 C12 — S23 Si3 S12 e"^ S23 Ci3 \ P , (2) 

\ iS iS I 

\ S23 S12 — C23 Si3 C12 e — S23 C12 — C23 Si3 S12 e C23 Ci3 y 

where P = diag(e~*'^^/^, e"*'^^/^, 1) contains the Majorana phases and Cij, Sij are defined as 
cos 9ij and sin6'jj, respectively. 

The smallness of neutrino masses is commonly attributed to the see-saw mechanism jS], 
which states that integrating out three heavy Majorana neutrinos with masses Mi^2,3 leads 
to a low energy neutrino mass matrix given by 

mj, = — m^M^^m^ . (3) 

Here M/j denotes the mass matrix of the three Majorana neutrinos A'^i 2,3 and m^i is a Dirac 
mass matrix resulting from the coupling of the Higgs doublet to the left-handed neutrinos 
and the Ni. 

On the other hand, the unexpected form of the mixing scheme implied by the data, namely 
one small (6'i3) and two large (^12,23 with 612 ~ ^'23) angles in combination with a rather 
moderate mass hierarchy, poses a more difficult problem. It might hint towards a particular 
symmetry, which seems not to be present in the quark sector. In order to find this under- 
lying symmetry of the neutrino mass matrix, one might consider minimal models and/or 
symmetries which are able to accommodate the data. For instance, one could introduce 
a non-standard fiavor charge corresponding to an Abelian ^7(1) symmetry and demand 
that the neutrino mass matrix conserves this fiavor charge. This will imprint a particular 
form onto rrii, and lead to certain mixing scenarios which can then be compared to the 
data. For three active flavor states there are the charges L^, and L^, out of which 
one can construct 10 linearly independent combinations of the form Cg Lf, + + Cr Lj. 
with Ce,^,T = or ±1. Neutrino data implies that these symmetries must be approximate. 
Nevertheless, already at the present stage we are capable of ruling out six of them, leaving 
us with four testable cases. Those four cases are easily to distinguish, because their pre- 
dictions for the neutrino mass spectrum and for the effective mass governing neutrinoless 
double beta decay |lj are very different: 

• conservation of implies the following form for the neutrino mass matrix j^]: 

/ \ 



m,. 



(4) 



\ - ■ dj 



2 



where a, b, d are in general complex parameters of order one. Introducing small break- 
ing terms in the first row of leads to a normal hierarchy, i.e., small neutrino masses 
rris ^ m2 > mi and a very small effective mass of order ^y^rn^ ~ 10~^ eV, where 



denotes the mass squared difference of solar neutrinos. The common mass 
scale of the matrix (jD) is the square root the mass squared difference of atmospheric 
neutrinos Am^; 

conservation of Lg — — demands the following form of the mass matrix jH] : 



/ a b\ 
■ 

V • ■ 0/ 



(5) 



The result is an inverted hierarchy, i.e., m2 ^ mi ^ m^ and an effective mass of order 
^y Am\ ~ 0.04 eV. Contributions from the charged lepton sector are crucial in order 
to reach accordance with the data because solar neutrino mixing is predicted to 
be maximal by Eq. Since the deviation from maximal neutrino mixing is of order 
of the Cabibbo angle , matrices conserving — — L^. are ideal candidates for 
incorporating Quark-Lepton-Complementarity (QLC) P; 

conservation of — has caught only little attention. The mass matrix implied 
by this symmetry has the form [TUl ITU IT^ : 



/ a \ 



m„ 



mo 



\ 



b 
■ 0/ 



(6) 



One is lead to quasi-degenerate light neutrinos, i.e., masses much larger than the 
mass splitting, and the effective mass is of order of the common mass scale mg ^ eV. 

Z^T, or Lf, Ln, Lf. Lj- or L^, + Ln Z^^, 



The other possibilities, i.e., conservation of L^, 
Lf. — + Lj- are ruled out phenomenologically. Finally, there is a tenth possibility cor- 
responding to the conservation of the total lepton charge + + L^-. The consequence 
would be that neutrinos are Dirac particles and that there is no neutrinoless double beta 
decay. Without more theoretical input, other than just + + Lr, the neutrino mass 
spectrum and the mixing angles are arbitrary (for some recent models for Dirac neutrinos, 
see for instance p!B1). 



As already mentioned, the probably most puzzling aspect of neutrino mixing is that one 
angle is close to zero and another one close to maximal. This indicates the presence of a 



3 



fi-T exchange symmetry ^1], which results in a mass matrix of the form 



f A B B\ 



rriy = 



D E 



(7) 



v • 



By comparing this matrix with the three matrices from Eqs. (jH EI) we see that only — 
Lr automatically possesses a /x-r symmetry. Therefore, 6*23 = 7r/4 and t/gs = are generic 
predictions of mass matrices obeying the flavor symmetry — Lr- Phenomenologically, 
tan ^23 is a priori not equal to one for (broken) and L^ — L^ — Lr-, but just a ratio of 
two numbers. 

On the more formal side, it turns out that the combination — is anomaly free and 
can be gauged in the Standard Model (SM) [Tbilisi- In fact, also the combination — 
or Le — Lr could be chosen to be gauged. However, (as indicated above) these two combi- 
nations are phenomenologically ruled out, since they typically predict two small neutrino 
mixing angles [T^ . 

Here we wish to analyze the consequences of a see-saw model based on the U{1) flavor 
symmetry — Lr- We will break the symmetry in a minimal way so that low energy 
neutrino data can be accommodated and that leptonic CP is possible. Interesting and 
testable phenomenology is predicted. Since neutrino physics has entered a precision era, a 
full study of any model has to take radiative corrections ^3 ^1 El 120] into account. In 
the framework of U{1) models under consideration, one should note the following: models 
based on Lg receive only minor corrections due to RG effects because these are in general 
very small for the normal hierarchy. If — — Lj. is the symmetry from which a 
model originates, then the opposite CP parities of the leading mass states suppress the 
running, thus leading again to small effects. However, for quasi-degenerate neutrinos, thus 
in particular for — L,-, one expects almost certainly large RG corrections. Moreover, a 
full study of RG effects in a see-saw model should not only consider the running from the 
lowest heavy neutrino mass Mi down to the weak scale, but also the running between the 
different heavy neutrino scales M3 to M2 and M2 to Mi [211 123 ESj • Those contributions 
can be at least as important as the usual ones from Mi to Mz- In our model we explicitely 
show that they are. 



Recently the case — Lr has been re-analyzed and a simple see-saw model has been 
constructed Apart from the Standard Model (SM) particle content there are only 

three extra SM singlets required. Those are heavy Majorana neutrinos which are used to 
generate light neutrino masses via the see-saw mechanism [H]. The charge assignment of 
the particles under L^ — Lr is given in Table [TJ As a consequence, the charged lepton mass 



2 See— saw Model for L 



L 



r 



4 













N2,m 


N3,TR 


$ 


L^i — Lj- 





1 


-1 





1 


-1 






Table 1: Particle content and charge under the U{1) symmetry corresponding to — L^-. 



matrix is real and diagonal and the Lagrangian responsible for the neutrino masses reads 

-C = {mD)ai + \NfC-^ {MR)ij Nj + h.c. 



$ {aVeNi + bV^N2 + dVrN^) + ^ {X NlC-^Ni + Y N^C-^N-i) + h.c. 



(8) 



Here, M is the high mass scale of the heavy singlets, $ the Higgs doublet with its vacuum 
expectation value v/ \/2 ~ 174 GeV and and C is the charge conjugation matrix. In terms 



of mass matrices, we have 





/ a 









f X 


\ 


mr, = V 





b 





and Mr = M 




Y 




V 





d) 




\ ■ 


■ / 



(9) 



At this stage, we are dealing with real parameters a, b, d, X, Y. After integrating out the 
heavy singlet states, the neutrino mass matrix = —rrijj is given by 



V 

'm 



( 

X 



\ 











\ 


bd 




Y 







/ 



(10) 



We need to break — in order to generate successful phenomenology^ . Note that the 
mass matrix from above has maximal 23 mixing, arbitrary 12 and zero 13 mixing as well as 
only one non-vanishing Am^, given by f ^/M^ [}? (ff' jY'^ — j X'^\ This is the mass splitting 
associated with the solar neutrino oscillations. We consider from now on the situation in 
which the symmetry is softly broken by additional small parameters in Mr. Moreover, in 
contrast to Ref. we allow for CP violation. Since the Dirac and charged lepton mass 
matrices are already diagonal, we can eliminate the possible phases of a, 6, d by redefining 
the lepton doublets. Consequently, prior to breaking only two phases are present, which 
we assign to the heavy singlets such that we replace in Eq. © X with X e*"^ and Y with 
Y e^^ . The first and minimal approach is to add just one small entry to Mr. For instance, 
we can add to the 12 element an entry ee*-^ with real e ^ 1. As a result, the low energy 



^For alternative breaking scenarios in case of CP conservation, see |12|. 
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mass matrix will read 



X 



V 

M 








XY 



\ 



bde-Z 
Y 

xV^ 



It is interesting to note that there is no CP violation in oscillation experiments or at high 
energy. This can be shown, for instance, by rephasing the elements of Eq. (fTT|) or by 
considering basis invariant quantities. Let us define here the quantity Jcp, to which any 
CP violating effect in neutrino oscillation experiments will be proportional. In terms of 
mixing angles and the Dirac phase 6 it is given by 

Jcp = lm{Uei U^2 = \ sin 2^i2 sin 2^23 sin 2^13 cos ^13 sin 5 . (12) 

From the mass matrix mjy, in turn, one can deduce [21] 

Im {hi2 /i23 ^31} = ^"^21 ^"^31 ^^32 Jcp ■, where h = mj, m,^ . (13) 

Using the mass matrix from Eq. (fTT|) to evaluate this expression yields that Jcp is zero. 
In order to have CP violation, we are therefore forced to add another perturbation to M^: 



'X e'' 



Mr = M 




(14) 



However, one can show that there is only one physical phase, and it therefore suffices to 
consider 

' X ei 

Mr = M \ ■ 626'^ Y \ (15) 



to obtain the following low energy effective mass matrix 

-adei 

X ^ 

,2 



rriu = - 



M 



XY 







hd 
Y 



\ 



-d' {X e2 - ej) 
XY' 



(16) 



For the rest of this contribution we will focus on this mass matrix. Let us therefore analyze 
in detail the implied phenomenology. 

A matrix with zero entries in the e/x and fifi elements has well-known predictions |25| |2S] ■ 
In particular, only quasi-degenerate light neutrinos are compatible with such a matrix 
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and in addition it is required that the ee and the /xr elements are leading and of similar 
magnitude [2^1 • This is however just the approximate form of a mass matrix conserving 

— Lr- 

By using the definition of the mass matrix from Eq. and inserting the conditions 
and m„„ = 0, one can obtain an expression for the ratio of the neutrino masses. 



m 



Expanding in terms of the small parameter lU^sl, one finds: 



mi 
ms 

m2 



tan^ 6*23 — lUesl cos 5 cot^^ 



tan 6- 



12 



23 



COS 623 



tan^ 6*23 + lUesl cos 5 tan 6*12 ^^^^23 

cos 



(17) 



23 



Since ^23 is close to maximal and ^13 is small, the light neutrinos are obviously quasi- 
degenerate. From the ratio of masses one can calculate the ratio of the solar mass squared 
difference over the atmospheric one: 



R 



AmQ ml — m\ 
Am 



^ \m^ — m. 



4 I Ue3 1 cos 6 



tan 6*23 sin^ 6*; 



cos 2^23 sin 26 



23 



'12 



As R is inverse proportional to the rather small quantity cos26'23, it is necessary that 
sin ^13 cos 5 = Ref/e3 <^ 1. Hence, the Dirac phase should be located around its maximal 
value 71/2, i.e., CP violation is close to maximal. The larger sin ^13 is, the smaller cos 5 
has to be, which indicates that for sizable \Ue3\ only maximal CP violation will occur. 
Furthermore, to keep R small, the angle ^23 can not become maximal. If it is located 
above 7r/4, then it follows from Eq. (fT7|) that mi > m^ and m2 > m^, which leads to an 
inverted ordering. Equivalently, for ^23 < 7r/4 one has m^ > mi and m^ > m2 and the 
mass ordering is normal. These aspects of the phenomenology are almost independent on 
the precise value of 612, which — as we will see later on — can receive large renormalization 
corrections. The effective mass governing neutrinoless double beta decay can be written 
as 

a'^ 9 ^ / Ami , , 

(m) = — - ~ tan2^23W^ , \^ , , (19) 



231 



which shows again that maximal atmospheric neutrino mixing is forbidden. For later use, 
we note here that for the neutrino texture under discussion typically ~ y92 — ti" holds, 
i.e., the two Majorana phases have very similar values [2n]- This will be very useful when 
we analyze the renormalization aspects of the model. 

Let us now take a closer look at the structure of the Dirac and heavy Majorana mass 
matrices. Since the charged leptons display a hierarchy, it is natural to assume that also 
the eigenvalues of the Dirac mass matrix are hierarchical. Then it is required that also the 
heavy Majorana neutrinos display a hierarchy in the form oiY ^ X. Typical values of 
the parameters which in this case successfully reproduce the neutrino data are Y = 0{1), 
a = C(O.Ol), 6 ~ ci = C(O.l), X = 0(0.001), ei = C(O.OOl) and 62 = O(O.l). With these 
values, the eigenvalues of Mr are approximatively given by MX and M {Y ± €2/2). We 
plot in Fig. ^ some of the correlations obtained for parameters around the values given 
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above. The oscillation parameters are required to lie within the So" ranges from Ref. P7] . 
In fact, just demanding zero entries in the e/x and /z/x elements will lead to such a behavior 
of the observables jSHj. Indeed, as can be seen from the Figure, the neutrinos are close in 
mass, atmospheric neutrino mixing is not maximal and lies above (below) 7r/4 for the in- 
verted (normal) mass ordering. Furthermore, CP violation is basically maximal for sizable 
values of \Uez\. 



3 Renormalization Group Effects 

As the flavor symmetry — L^. leads to quasi-degenerate neutrino masses, strong running 
of the mixing angles is generically expected [18,, .20j. Typically, the running of the mixing 
angles 9ij in a quasi-degenerate mass scheme with a common mass scale mo is proportional 
to m1/{mf — rrij) and therefore particularly strong for 612- This behavior turns out to hold 
not only for the running between the scale of the lightest heavy Majorana neutrino and 
the weak scale, but also when the running between the see-saw scales is taken into account 
PT| 1^ . The heavy Majorana neutrinos have to be integrated out one after the other, 
leading to a series of effective theories. It turns out that the running of the observables 
between the scales M3 and M2 or M2 and Mi can be of the same order as the running 
between Mi and Mz- In general this leads to quite involved expressions for the /3-functions. 
However, in our case the structure of the Dirac and charged lepton mass matrices (i.e., the 
fact that they are diagonal) simplifies matters considerably and allows for some analytic 
understanding of the numerical results. All of the expressions we give in this Section are 
understood to be precise to order \Ue3\ and can be obtained with the help of Refs. [2II1I22]- 
For instance, the /3-functions for the mixing angles are: 

2 

167r2 9u (1 + cos ((^2 - <^i)) sin 2^12 [Pu - {P22 cos^ 623 + P33 sin^ ^^23)] , 



327r2^^i3--r^(cos((5-v^i)-cos(5-^2))sin2^i2sin2^23 (P22 - P33) , (20) 

2 

167r^e23 [(1 + cos(/?2) cos2^i2 + (1 + cos(/?i) sin^^ia] sin 2^23 (^22 - ^^33) • 

/\Tn,j^ 

Note that since ipi ~ ip2 there is no cancellation in the first relation for ^12. The matrix P 
in the above equations is given in the SM and in the Minimal Supersymmetric Standard 
Model (MSSM) by 

[ Y^Y} + r^rj MSSM . 

Since in our model both Fj, = mo/v and Yf. = mg/v are diagonal, P is diagonal as well 
and no off-diagonal elements appear in Eq. (f^UI) . With the values of the parameters in 
rriD as given in Section |21 (i.e., a <^ h ^ d = (9(0.1)) we can safely neglect Pu. This 



o 

2 
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leads in particular to a negative /^-function for 612 ■ Hence, this angle will become larger 
when evolved from high to low scales. Neglecting further the charged lepton Yukawas in 
Ye above the see-saw scales and noting that P22 — &^/2 and P33 ~ for the SM and 
twice those values for the MSSM, we see that the running of 6'i3 and 623 is suppressed with 
respect to the running of 612 due to two reasons: firstly, it is inverse proportional to Am,\ 
and secondly, it is proportional to P22 ~ -P33 cc b"^ — (P, which is smaller than P22 + P33, 
on which the running of 612 (approximately) depends. Hence, the running of 6*13 and 6'23 is 
suppressed by R (6^ — d?)/b'^ above the see-saw scales. 

Below the see-saw scales only the tau-lepton Yukawa coupling y^- governs the RG correc- 
tions. In this regime the evolution goes like pUj 

2 

16 TT^ 9^2 ^ -vl sin 29^2 sin^ (1 + cos(yP2 - Vi)) , (22) 

in the MSSM, which is again negative and leads to sizable running. The formulae for the 
running of 6*13,23 are again suppressed by roughly a factor R. 

The phases stay almost constant in the whole range, because one can show that for our 
model ipi and 02 are mainly proportional to sinyji and siny92, respectively. Roughly the 
same behavior is found below Mi. Finally, the RG effects on the neutrino masses correspond 
predominantly to a rescaling, since the flavor-diagonal couplings, i.e., gauge couplings and 
the quartic Higgs coupling, dominate the evolution [T^ . 

We can analyze if the zero entries in the mass matrix Eq. remain zero entries. Below 
the see-saw scale it is well known that the RG corrections are multiplicative on the mass 
matrix, a fact which leaves zero entries zero. Taking the running in between the heavy 
Majorana masses into account, one can factorize the renormalization group effects Zext 
from the tree-level neutrino mass matrix in the MSSM as 



= Zj^^m^Z^^t . (23) 

Since in our model the neutrino and the charged lepton Yukawa coupling matrices are 
diagonal, the renormalization group effects are flavour diagonal, too. Therefore, texture 
zeroes in the charged lepton basis remain zero. With the already mentioned simplifications, 
Zext is approximately^ given by 



,2, MY MY 

+ b^ In — h (T In — h 

A A 



- 3^2 + 3z/i 



In^, (24) 



where ((^2, 5^1) are the gauge couplings^ of the electroweak interactions SU{2)i x f/(l)y, y^ 
and yr are the top quark and r Yukawa couplings, respectively and A is the GUT scale of 
2 ■ lO^^ GeV. 



^As the perturbations are small, the mass eigenvalues of the heavy right-handed neutrinos are approx- 
imately given by (Afi, M2, Afg) ^ M {X, Y{1 - 62/2), y(l + £3/2)) ~ {X, F, Y). 

■^We use GUT charge normalization for gi, i.e., gf^ — y fffi- 
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L' 



matrix 



comments 



normal hierarchy 
always small RG effects 




R 



AniQ 



\U, 



e3\ 



tan2 6'atm ^l + \Ue3 

(m) : 



v/Ami|f/,3|2 



1 + Vr 



m 



Q 



Le ~ — Lt 

inverted hierarchy 
small RG effects (pseudo-Dirac) 




requires Ui>: ideal for QLC 

tan^^Q ~ 1 - 4|f/e3| ~ 1 - 2^2 sin 
(m) ~ >y Am\ 



— Lj- 
quasi-degeneracy 
always large RG effects 




in leading order: 
Ue3 = and 6*23 = 7r/4 
(m) ~ tuq 



Table 2: The three simple U{1) charges as implied by current data, the zeroth order mass 
matrix and some comments on their phenomenology. 



In the extended SM, however, there are additional corrections which can not be factorized. 
They are responsible for the instability of texture zeroes under the renormalization group 
in the SM |^ I23j. We checked that for most observables the running behavior in the SM 
is similar to the running in case of the MSSM. The solar neutrino mixing angle receives 
however more RG corrections in the SM, a fact which can be traced back to the appearance 
of an e/i entry in the mass matrix Eq. ()16|) . One might wonder at this point if this filling 
of zero entries would allow one to generate successful phenomenology from a mass matrix 
obeying the flavor symmetry — Lj. without any breaking, i.e., just from Eq. (jU)). Recall, 
however, that in the SM — is anomaly free and therefore the texture of the mass 
matrix is stable. 

We plot in Fig. |21 the running of the angles, phases and masses for a typical example in 
the MSSM for tan (3 = 10 [28 . The resulting neutrino parameters at the GUT scale are 
sin^^ia = 0.123, Ue3 = 0.0484 e*^•^^ sin^ ^23 = 0.481, (mi, m2, ms) = (0.1527, 0.1533, 0.1653) eV 
with AttIq = 1.9 ■ 10~^ eV^ and Am\ = 3.8 ■ 10~^ eV^. They are changed by the renor- 
malization group evolution to sin^^i2 = 0.303, U^s = 0.0496 e^^■^^ sin^ ^23 = 0.481, 
(mi, m2, m3) = (0.1152, 0.1155, 0.1245) eV with Am| = 7.9 • 10"^ eV^ and Am^ = 
2.1 ■ 10~^ eV^. Note that the phases and ^13,23 remain practically constant, whereas sin^ 6'i2 
and Am^ ^ are changed by factors of up to three, and that the running in between and 
above the see-saw scales is at least as important as the running below them. This implies 
that radiative corrections, in particular for 6'i2 and the mass squared differences, can be 
crucial if one wishes to make meaningful phenomenological predictions of a given model. 
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4 Conclusions and Summary 

Identifying the structure of the mass matrix m,y is one of the main goals of neutrino physics. 
Current data seems to point towards three interesting forms of rriy, each corresponding to 
the (approximative) conservation of a non-standard lepton charge. While and Lg — L^ — 
Lj. are well known cases — leading to the normal and inverted hierarchy, respectively — the 
third possibility — Lr has not been examined in as much detail. Since quasi-degenerate 
neutrinos are predicted by the conservation of L^—Lr, this case is easily testable. Moreover, 
since the mass matrix is automatically //-r symmetric, further phenomenological support 
is provided. We stress again that the three possibilities Lg, L^. — L^ — and — Lj. differ 
in particular by their prediction for neutrinoless double beta decay, a fact which will easily 
allow to identify the correct one. In Table|21the relevant matrices and some aspects of their 
phenomenology are summarized. Here we analyzed the consequences of a see-saw model 
based on — Lr, which predicts, e.g., non-maximal atmospheric neutrino mixing, large 
CP violation and quasi-degenerate neutrinos. This last aspect makes radiative corrections 
be particularly important and we carefully analyzed their effects both numerically and 
analytically. We showed explicitely that running between the see-saw scales is at least as 
important as running from the lowest see-saw scale down to low energies. 
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Figure 1: Scatter plots of some of the observables resulting from a mass matrix with zero 
e/z and fifi elements. We required the 3a values for the oscillation parameters from |i27j . 
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Figure 2: Typical plot for the evolution of the mixing angles and masses from the GUT scale 
2 X 10^^ GeV to low scale. In this example, we have tan (3 = 10 and the parameters in the 
matrix Eq. ^ at the GUT scale are a = 0.0066926, b = 0.0692883, c = 0.0697464, X = 
0.0096528, y = 1, ei = 0.0005595, es = 0.0749098, ip = 2.45376 and M = 9.098937108-10" 
GeV. 
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